Dirac Matrices for Chern— Simons Gravity 
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A genuine gauge theory for the Poincare, de Sitter or anti-de Sitter algebras can be constructed 
in (2n — l)-dimensional spacetime by means of the Chern-Simons form, yielding a gravitational 
theory that differs from General Relativity but shares many of its properties, such as second order 
field equations for the metric. The particular form of the Lagrangian is determined by a rank n, 
symmetric tensor invariant under the relevant algebra. In practice, the calculation of this invariant 
tensor can be reduced to the computation of the trace of the symmetrized product of n Dirac Gamma 
matrices T a t in 2n-dimensional spacetime. While straightforward in principle, this calculation can 
become extremely cumbersome in practice. For large enough n, existing computer algebra packages 
take an inordinate long time to produce the answer or plainly fail having used up all available 
memory. In this talk we show that the general formula for the trace of the symmetrized product 
of 2n Gamma matrices T a b can be written as a certain sum over the integer partitions s of n, with 
every term being multiplied by a numerical coefficient a 3 . We then give a general algorithm that 
computes the a-coefficients as the solution of a linear system of equations generated by evaluating 
the general formula for different sets of tensors B ab with random numerical entries. A recurrence 
relation between different coefficients is shown to hold and is used in a second, "minimal" algorithm 
to greatly speed up the computations. Runtime of the minimal algorithm stays below 1 min on a 
typical desktop computer for up to n = 25, which easily covers all foreseeable applications of the 
trace formula. 



I. INTRODUCTION 

There's more to higher-dimensional gravity than Ein- 
stein and Hilbert [J-Q. 

Chern-Simons (CS) gravity in d = 2n— 1 dimensions is 
a gauge theory for the Poincare, de Sitter or anti-de Sitter 
(AdS) algebras, depending on the value of the cosmolog- 
ical constant [1, [^] . 

Let us focus on the AdS algebra, SO (d — 1,2). A 
convenient matrix representation is provided by Tab = 
F^Ts], where Ta are Dirac matrices in D = d + 1 = 2n 
dimensions 



where {• • • } denotes symmetrized matrix product. 
The pseudoscalar trace reads 

Tr(r„, {T Al B! ■ ■ •r.4„B„}) = 7 e AiBi— A n B n , 



(3) 



where 7 is a numerical coefficient. Use of this invariant 
polynomial brings in the Lanczos-Lovelock 0, Q family 
of Lagrangians into CS gravity. 

The scalar trace, on the other hand, is more involved. 

In this work we provide two algorithms that can be 
used to efficiently compute the scalar trace for any n 
and for any spacetime dimension d (without any implied 
relation between n and d). 
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(i) 

(2) 



The Lagrangian for CS gravity is shaped to a great 
extent by a rank-n, AdS-invariant symmetric polynomial 

(•••)■ 

This polynomial can be identified with any of the fol- 
lowing traces: 

• Tv {Ta 1 b 1 • • '^A n B n } (Lorentz scalar) 

• Tr (r» {TAtBr ■ ■ • ^A n B n }) (Lorentz pseudoscalar), 
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II. FORMULATION OF THE PROBLEM AND 
RESULTS 

Let us consider Dirac matrices T a , a = 0, . . . , d — 1, in 
rf-dimensional Minkowski spacetime. By definition, they 
satisfy the Clifford algebra [l(| 



T n Th + TfTn 



2r) ab l, 



(4) 



where r\ a \, = ( — !-•••+) is the usual Minkowski metric 
and 1 stands for the mxm unit matrix, with m = 2 L d / 2 J . 

The T-matrices which are the subject of this work are 
defined as 



r a b — r[ a r fc ] — - (r a r b - r b r a ) . 



(5) 



For completeness, let us define the symmetrized prod- 
uct of n matrices Mi, i = 1, . . . , n, as 

{M x ■ ■ ■ M n } = — • • • M An) , (6) 
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where the sum extends over all permutations tt in the 
symmetric group S n . 

Experience shows that the trace is most efficiently writ- 
ten with all matrices multiplied by arbitrary antisym- 
metric tensors. Take, for instance, the trace of the sym- 
metrized product of two Gamma matrices, and compare 
the following equations: 

Tr {r ah r cd } = m (rjadVbc - VacVbd) , (7) 
A ab B cd TT {TabTcd} = 2mA a b B b a . (8) 

The two terms on the right-hand side of eq. ([7]) have 
collapsed into one in eq. ([5]). Greater simplifications are 
achieved for more complicated cases. If desired, eq. ([7]) 
can be recovered from eq. © by means of the formal 
replacement A ab ->• 5f d , B ab ->• 8f d , where 5f d is the 
generalized Kronecker delta. 

Let Bf , i = 1,2,3,..., be arbitrary antisymmetric 
tensors, and let us define 

ft = BfT ab . (9) 

The symmetrized product of n /3-matrices can be 
written as a linear combination of matrices r ai ... = 
r [oi-" r ap]) witn P = 0, 4, 8,..., 2n (for n even) or 
p = 2, 6, 10, . . . ,2n (for n odd). The only term that con- 
tributes to the trace is that proportional to the identity 
matrix (p = 0). For odd d, however, the r 0l ... 0d matrix 
is also proportional to the identity and must be generi- 
cally taken into account when computing the trace. The 
expansion of the symmetrized product of the /3-matrices 
includes only T-matrices with an even number of indices, 
so that the r ai ... a<J -term never actually shows up in our 
case. In particular, this means that the trace of the sym- 
metrized product of an odd number of /3-matrices van- 
ishes identically. 

The trace of the symmetrized product of 2n /3- 
matrices, on the other hand, can be written as 

= m£a.B« (10) 

where the notation shn [Tlj indicates that the sum must 
be performed over all integer partitions s of n, and 
stands for the following sum of contractions of B-tensors: 

B«= £ f[(B i2 , i+ ... +2aj _ i+1 ---B i2ai+ ,.. +2Sj ). 

(il— %an) J =1 

.(11) 

In eq. (|11[) . the notation (i\ ■ ■ 'ian) is used to indicate 
that the sum must be performed over all «i,...,«2n £ 
{1, . . . , 2n}, with the restriction that they be all differ- 
ent. This implements the permutation of all /3-matrices. 
Every term in the sum contains the product of r factors 
of the form (Bi ■ ■ ■ B q ), where r is the length of the par- 
tition s = (si, . . . , s r ), with n — s\ + ■ • • + s r . The j-th 
factor in the product represents the trace of the product 
of 2s j S-tensors, i.e., 

(B 1 ---B q ) = (B 1 ) c \ 2 (B 2 r c3 .--(B q r ci , (12) 



with q = 2s j. 

To every term in eq. (jlOp . i.e., to every partition s of 
n, there corresponds an a s coefficient. Numerical values 
for the a-coefficients corresponding to the partitions of 
n = 1 , . . . , 7 are given in Table |TJ 

The following examples for n = 1 , . . . , 4 should help 
clarify the meaning of eqs. (fTU)) and (fTTj) : 

Tr {0102} = m^2a 1 (B i B j ), (13) 

Tr {0 1 ---/3 4 } = mJ2 ["2 (BiBjB k Bi) + 

(ijkl) 

+a u (BiBj) (B k Bt)] , (14) 

Tr{0 1 ---0 6 } = m [as {B X1 ■ ■ ■ B l6 ) + 

+ a 2 i (B h ■ ■ ■ B^) (B i5 B i6 ) + 

+am {Bi t B l2 ) {B l3 B u } {B l5 B l6 }] , (15) 

Tr{/3 1 ---/3 8 } = m ^ [a 4 (B h ■ ■ ■ B ia ) + 

+ a 3 i (B i± ■ ■ ■ B i6 ) (B i? B is ) + 

+ a 22 (B n ■ ■ ■ B u ) {B l& ■ ■ ■ B l8 ) + 

+ ol 2 \\ (B il ■ ■ ■ Bi 4 ) (Bi s Bi 6 ) (Bi 7 B ia ) + 

+ami {B n B i2 } (B i3 B iA ) (B is B i6 ) (B l7 B is )] . 

(16) 

The proof of eq. (|10[) is by exhaustion; the right-hand 
side includes all possible terms that may contribute to 
the trace of the symmetrized product of 2n /3-matricesQ 

Our approach to the computation of the a-coefficients 
is the subject of section HTT1 

III. METHOD 
A. General Algorithm 

The central observation behind the algorithm used in 
the computation of the a-coefficients shown in Table U is 
the fact that eq. (TIT))) is valid for arbitrary tensors Bf . 

For illustration purposes, let us focus first on the n = 3 
case. Eq. (fT5"j) simplifies greatly if we choose all B-tensors 
to be equal, since in this case the sum over all different 



2 The formula for Tr (T, {/3i • ■ ■ /3n}) includes pseudoscalar terms 
that appear in certain dimensions d (e.g., t a bcd^i b Bj d for d = 4) 
but are absent from Tr{/h • • • /?2n}, where only Lorentz scalars 
are allowed. Here T* = Tq ■ ■ ■ T^—i is the d-dimensional gener- 
alization of 75 in d = 4. 
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TABLE 


I. ^-coefficients corresponding 


to the partitions of 


n = 1, . . 


.,7. 




n 


s 


a s 


1 


1 


1 


2 


1 + 1 


1/2 




2 


—2/3 


3 


1 + 1 + 1 


1/6 




2 + 1 


-2/3 




3 


32/45 


4 


1+1+1+1 


1/24 




2 + 1 + 1 


-1/3 




2 + 2 


2/9 




3+1 


32/45 




4 


-272/315 


5 


1+1+1+1+1 


1/120 




2+1+1+1 


-1/9 




2 + 2 + 1 


2/9 




3 + 1 + 1 


16/45 




3 + 2 


-64/135 




4+1 


-272/315 




5 


15872/14175 


6 


1+1+1+1+1+1 


1/720 




2+1+1+1+1 


-1/36 




2+2+1+1 


1/9 




2 + 2 + 2 


-4/81 




3+1+1+1 


16/135 




3 + 2 + 1 


-64/135 




3 + 3 


512/2025 




4+1 + 1 


-136/315 




4 + 2 


544/945 




5 + 1 


15872/14175 




6 


-707584/467775 


7 


1+1+1+1+1+1+1 


1/5040 




2+1+1+1+1+1 


-1/180 




2+2+1+1+1 


1/27 




2+2+2+1 


-4/81 




3+1+1+1+1 


4/135 




3+2+1+1 


-32/135 




3 + 2 + 2 


64/405 




3 + 3 + 1 


512/2025 




4+1+1+1 


-136/945 




4 + 2 + 1 


544/945 




4 + 3 


—8704/14175 




5 + 1 + 1 


7936/14175 




5 + 2 


-31744/42525 




6 + 1 


-707584/467775 




7 


89473024/42567525 



permutations of i±, . . . , i§ £ {1, . . . , 6} is trivially per- 
formed. The result reads 

J-Tr (^) = a 3 (-B 6 ) + a 21 (S 4 ) (S 2 ) + a ni (B*f . 

(17) 

We wish to cast eq. (fl7|) as a linear equation with three 
unknowns, namely, a 3 , a 2 \ and am- To do this we need 
to be able to assign numerical values to the left-hand 
side and to the various (B q )-terms that appear on the 
right-hand side. We accomplish this by (i) picking some 
antisymmetric tensor B ab with random numerical entries 
and (ii) choosing an explicit representation for the T- 
matricesH We emphasize that the possibility of choosing 
the i3-tensors at will relies upon the fact that eq. (fTOj) is 
valid for arbitrary B^s. 

To be able to solve for the a-coefficients we need two 
more equations. These are readily obtained by randomly 
selecting two further B-tensors. Denoting the three dif- 
ferent choices for the B-tensors by B/,, with k = 1,2,3, 
we obtain the following 3x3 linear system: 



-7(111) 

Z{ 'am - 


h Z\ 'a 21 - 


1 -7(3) 
r- Z{ 'a 3 


= T U 


(18) 


,7(111) 

Z\ 'am - 


h Z\ 'a 2 i - 


h zf ] a 3 


= T 2 , 


(19) 


^3 Cfclll " 


v( 21 ) 
V Z\ 'a 21 - 




= n, 


(20) 



where 

Tk = ^-Tr (0°) , (21) 
6!m 

Z^ = (Bl) 3 , (22) 

= (Bt) (Bl) , (23) 

Zf = (Bt) • (24) 

The method to compute the a-coefficients for any value 
of n is now clear and can be summarized in the following 
sequence: 

1. Let p = p (n) be the number of partitions of n0 

2. Choose an explicit representation for the T- 
matrices (see, e.g., Ref. 

3. For k — l,.,.,p, do: 

(a) Pick an antisymmetric tensor B% b with ran- 
dom numerical entries. 

(b) Compute 

T k = -\ r Tv{Pl n ) 1 (25) 
(Zn)\m 

where k = BfT ab . 



3 See section ITVl for a discussion of the choice of spacetime dimen- 
sion d in which to carry out the computation. 

4 The function p (n) is called the "partition function" in the math- 
ematical literature [Til . 
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(c) For every partition shn, with n = si 
s r , compute 



As) 



n K j 



(26) 



The notation (i?^) stands for [see eq. (fl2)) ] 
(^) = (5 fe ) Cl C2 (B fe ) C2 C3 ---(B fe ) c " ci . (27) 

4. The a-coefficients are the solution to the pxp linear 
system of equations 

X Z /~"- ll (k = l,...,p). (28) 



B. Minimal Algorithm 

Careful inspection of the a-coefficients shown in Ta- 
ble U shows that there exists a recurrence relation among 
different coefficients. 

Let s be a partition of n. The frequency representa- 

), where fij 



tion 



ll| of s is the notation s = (l" 1 2^ 2 
represents the multiplicity of j, i.e., the number of times 
that a given integer j appears in s. 

We find that the coefficient a s corresponding to the 
partition s — (l Ml 2^ 2 • • • ) can be written as 



J'j 



(29) 



where ctj are the coefficients associated with the "elemen- 
tary" partitions f = f , 2 = 2, 3 = 3, etc. 

For example, all coefficients associated with the non- 
elementary partitions of n = I , 2, 3 can be computed from 
ot\, a.2 and a.3 by means of the equations 



an 


a\ 
21 


«§ 

0! 




1 

2' 




(30) 


" 1 ' 1 = 


¥ 


«§ 
0! 


0! 




1 

6' 


(31) 


C*21 = 


IT 


"5 
1! 


«§ 
0! 




2 

_ 3' 


(32) 



Of course, this recurrence relation also holds for more 
complicated cases, such as 



Q!3211 



a\ a\ a\ a\ a\ a% a° 

2T U U "oT"oT"oT"oT 



32 
135 



(33) 

29]) y ields 



When applied to an elementary coefficient, eq. 
an identity. 

The recurrence relation in eq. (I2U1) can be used to com- 
pute the values for the a-coefficients associated with all 
the non-elementary partitions of n. Its use, however, re- 
quires knowledge of the elementary coefficients, for which 
no closed formula is available. This situation suggests a 



"minimal" algorithm that (i) calculates elementary coef- 
ficients in a manner analogous to that of the "general" 
algorithm and (ii) computes non-elementary coefficients 
from eq. (j2"9"j) . 

The following sequence describes such an algorithm: 

1. Let N be the maximum integer for which we wish 
to calculate the a-coefficients. 

2. Choose an explicit representation for the T- 
matrices. 

3. Pick an antisymmetric tensor B ah with random nu- 
merical entries 

4. For n — 1, . . . , N, do: 

(a) Compute 

1 



T = 



(2n)\m 



Tr(/3 2 



(34) 



where f3 = B ab T ab 



(b) For every partition shn, with n = si 
s r , compute 

= f[(B 2s >). 

i=i 



(35) 



(c) Use the recurrence relation (f29|) to calculate 
all non-elementary coefficients associated with 
the partitions of n (this step is empty for n = 
!)• 

(d) Solve 



Y J Z {s) a s =T 



(36) 



for a n (this is a linear equation with one un- 
known) . 



IV. DISCUSSION AND CONCLUSIONS 

The algorithms described in section IIIII turn around 
the problem of finding formulas for the trace of a prod- 
uct of Gamma matrices. The usual textbook approach 
starts with eq. @ and deduces the required formulas 
from there. Our approach here works the other way 
around. We start by identifying the general form of the 
equation for the trace of the symmetrized product of 2n 
/3-matrices. Eq. (|10p amounts to such an identification, 



We took d = 2 and B 01 = +1, since a two-index antisymmet- 
ric tensor has only one degree of freedom in two spacetime di- 
mensions, and overall numerical factors are not significant for 
the calculation. See section IIVI for a discussion of the choice of 
spacetime dimension d in which to carry out the computation. 
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since it contains all possible sums of ^-contractions that 
may contribute to the trace. The a-coefficients appear 
as undetermined parameters, which are computed by de- 
manding validity of eq. (|10[) in several nontrivial cases. 

As stressed in section [TTTl our method works because 
cq. (|10p holds for arbitrary antisymmetric tensors Bf b . 
We have used B-tensors with random numerical entries 
to generate the linear system of equations whose solu- 
tion provides the a-coefficients. In this sense our ap- 
proach bears some resemblance to Monte Carlo methods, 
where random numbers play a crucial role. The use of 
random matrices^ however, is not essential to our cal- 
culation. All that is required for the general algorithm 
to succeed is a set of -B-tensors such that every iteration 
produces an equation for the a-coefficients that is lin- 
early independent from the rest, yielding a full-rank Z 
matrix [cf. eq. ([28]) ]. 

The solution we find is, of course, independent of the 
choice of B-tensors; this is conceptually clear, but can 
also be verified by running the algorithm several times 
with different sets of (randomly generated) B-tensors. 
The fact that the same solution is obtained every time 
confirms both this independence and the correctness of 
eq. (flO]) . i.e., that no other terms can be added to the 
trace. 

The a-coefficients are also independent of the space- 
time dimension d, which means that the algorithm should 
in principle work for any d we choose. There is, however, 
an important caveat. To produce a solvable system one 
needs the _B-tensors to have a sufficient number of inde- 
pendent components, so that the successive iterations of 
the algorithm yield linearly independent equations. We 
find that there is a minimum spacetime dimension d = 2n 
that allows the Z matrix to achieve full rank. This means 
that the general algorithm must be run with d > In in 
order for a solution to be produced. 

The minimal algorithm, with only one linear equation 
to be solved, works even with a minimum spacetime di- 
mension of d = 2. 

Is our approach any better than the textbook method? 
One way to probe into this question is to compare the 
runtime of both. The textbook method can be imple- 
mented in, e.g., Kasper Peeters' excellent computer alge- 
bra system "Cadabra" [l3l Hij, We were able to deduce, 
starting only from the definition of Dirac matrices, the 
a-coefficients for n = 1,2,3. The n = 3 case took some 
30 min to be solved on a typical desktop computer^ while 
the n = 4 case caused the program to crash. This ap- 
proach, of course, requires hardly any input and produces 
the full sought-after formula. Starting from eq. (fTU|). we 
programmed our general algorithm in the computer al- 
gebra system "Maxima" [l5f and were able to run it suc- 



cessfully for n = 1, . . . , 7. The n = 8 case caused Maxima 
to run out of memory. Runtime for n — 1, . . . , 4 was neg- 
ligible, while the n = 7 case took under half an hour. The 
minimal algorithm, which we also programmed in Max- 
ima, had negligible runtime even for N = 25. Table |TT] 
summarizes runtime for these different scenarios. 

Complexity for the general algorithm grows exponen- 
tially with n. Complexity for the minimal algorithm, on 
the other hand, grows linearly with p0 All foreseeable 
applications of the formula for the trace of a product 
of In Gamma matrices are well covered by the minimal 
algorithm with negligible runtime. 

TABLE II. Approximate runtime for the textbook method 
(as implemented in Cadabra) and the general and minimal 
algorithms (as implemented in Maxima) on a typical desk- 
top computer. For the minimal algorithm, the first column 
is understood to mean N, the maximum integer for which 
the a-coefficients are computed. The second column lists the 
partition function of n, which corresponds to the number of 
a-coefficients to be determined. 



n 


V 


Textbook 
Method 


General 
Algorithm 


Minimal 
Algorithm 


1 


1 


negligible 


negligible 


negligible 


2 


2 


negligible 


negligible 


negligible 


3 


3 


~ 30 min 


negligible 


negligible 


4 


5 


crashed 


negligible 


negligible 


5 


7 




negligible 


negligible 


6 


11 




~ 1 min 


negligible 


7 


15 




~ 24 min 


negligible 


8 


22 




crashed 


negligible 


9 


30 






negligible 


26 


2436 






~ 1 min 


28 


3718 






~ 2 min 


30 


5604 






~ 5 min 
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6 To be precise, what we use are actually two-index antisymmetric p (n) ^ given by thg H ardy-Ramanujan equation, p (n) 
tensors with random numerical entries. 

7 In 2011 this meant a 3.20-GHz CPU, with 3.7 GB of memory. 

8 An asymptotic approximation for the partition function 



l/4n v / 3) cxp (n^/2n/3) 
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